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Abstract 



The branching ratio for the ;U — )• 67 decay in the framework of the minimal 
flavour violation in the MSSM is calculated for various regions of the MSSM pa- 
0^, rameter space. The lepton flavour violation goes through the PMNS mixing matrix. 

p, ^1 The dependence on tan (3 is studied in comparison with experimental data. The re- 

O ■ suits crucially depend on the mixing angle ^13. Observation of this decay would 

serve as a manifestation of new physics beyond the SM. 



: 1 Introduction 
00 , 

The lepton flavour violating (LFV) processes are forbidden in the Standard Model if all 
^ . neutrinos are massless. They neither go via the charged currents nor via neutral ones. The 
Q ! situation changes if neutrinos are massive. Then one has to introduce the right-handed 
neutrinos to the SM and the lepton sector becomes similar to the quark one. By analogy 
with the CKM mixing matrix one obtains the PMNS mixing matrix whose elements are 
measured today with a reasonable accuracy. The LFV processes become possible via the 
^ . charge currents. FCNC are still forbidden. 

, The /i — 67 decay is one of the most important LFV processes. In the SM with 

neutrino masses it is described by the so-called penguin diagram and is proportional to 
the off-diagonal elements of the PMNS matrix. Due to chirality conservation it is also 
proportional to the neutrino masses. Calculations in the SM give [lH5] 



BrsM (/i 67) 



3a 
32^ 



With the present neutrino oscillation data one can estimate the branching ratio BrsuifJ' 
67) ~ 10~^^, a ridiculously small number. So the only reasonable possibility to observe 
such a decay would be a non-minimal flavour violation or a new physics with heavy 
particles propagating in the loops. Modern experiments [B] have reached the upper bound 
of 

Sr(/i^e7) < 2.8- 10"^\ 
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and the MEGA experiment [7j gives a similar result 

Br{fi 67) < 1.2 ■ 10'^\ 

They observe some "suspicious" events, though of no statistical value so far. 

One sees that the experimental bound is far from the SM prediction. However, in ex- 
tensions of the SM, for instance, in the Minimal Supersymmetric Standard Model (MSSM) 
there are additional contributions to the LFV processes. In particular, one can easily get 
the branching ratio for the fi ej decay in the MSSM of the order of 10~^^, which makes 
the measurement of this decay rate extremely sensitive to new physics. One can say that 
if the decay is observed this would be the manifestation of physics beyond the SM. 



2 Mixing matrix in the lepton sector 

When neutrinos are massive, the lepton sector of the SM resembles the quark one and 
naturally includes the mixing of flavours. The left-handed fields of the neutrino flavour 
eigenstates i/e^, ^fi^, s-re linear combinations of the three neutrino mass eigenstates 

3 

i=l 

and the weak interaction via the charged lepton current reads 

Ccc = E 'iLlaiyi.W''^ + h. c, (2) 

where 



U= \ f/^i U^2 (3) 





Ue2 


Ue3 


U,l 


U,2 


U^3 


Url 


Ur2 


Ur3 



is the unitary neutrino mixing matrix, the so-called PMNS matrix [8l[9]. The PMNS 
matrix can be parametrised by three angles and, depending on whether the massive 
neutrinos are Dirac or Majorana particles, by one or three CP violating phases [TUHT^ . 
We write 

(C12C13 S12C13 sue'^ \ 

^■5l2C23 — Ci2Si3S23e*'' C12C23 — ^123135236*'^ C13S23 T (4) 
iS iS I 

S12S23 ~ Ci2Si3C23e — C12S23 — C23Si2Si3e C13C23 / 

where Cij = cos 6ij, Sij = sin 6ij, 6 is the Dirac CP violating phase, and the matrix 



T = di&g[l,e'^,e'^j (5) 

contains the Majorana CP violating phases 021 and 0:31. In the following, we neglect all CP 
violating phases assuming CP conservation in the lepton sector. The existing experimental 
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neutrino oscillation data allow the determination of the solar and atmospheric neutrino 
oscillation parameters 612 and 623 with a relatively good precision. It is also possible to 
place rather stringent bounds on the angle ^13. By analysing the experimental data it has 
been found that [T8|[Ti] 

sin^^i2 = 0.304in?6, sin2^23 = 0.5^°;°^. 

Even though the angles 612 and 623 are known with a reasonable accuracy, the angle 
613 still remains unknown and is the main source of uncertainties of our predictions. A 
combined 3-neutrino oscillation analysis of the global data gives [H] 

sin^^is < 0.035 (0.056) at 90% (99.73%) C.L. 

and a global analysis of all available neutrino oscillation data provides the numerical 
value [15] 

sin^ = 0.016 ±0.010. 

As it will be clear later, the — )■ 67 decay rate is proportional to sin6'i3 and vanishes 
with the latter. 

3 The /i — > 67 decay rate 

The main decay of muon is the double neutrino decay /i eu^Ue which gives almost 
100% of the width and is one of the best measured decays. In the SM it goes through the 
W^-boson exchange (see FigH]) and the decay width is given by 




Figure 1: The double neutrino muon decay 



where 

F{x) = 1 - 8x - I2x^ In X + - (7) 

accumulates the radiative corrections and the Fermi constant is related to the weak cou- 
pling and the PF-mass by Gp = g'^ / 4:\/2m^ . Since rrie ^ m^, the value of F(ml/m'^) = 
0.9998 is very close to 1. 

The vertices in FigJT] contain the PMNS matrix elements. In the case when one does 
not distinguish between neutrino flavours, one has to sum over all neutrino species and 
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due to absence of interference obtains the sums J^i^f^i^U ~ ^ Ylii^ei^ei = 1- The 
latter are a consequence of the unitarity of the PMNS matrix. So the decay width ^ 
coincides with the one obtained in the SM without neutrino mixing. On the contrary, 
if one distinguishes neutrino flavours, the decay rate is proportional to specific PMNS 
matrix elements and is smaller. 

Jl!h/\an Jr^n^ JLtM^ 



IJ V e ^ V e y ^ 

a) b) c) 

Figure 2: The penguin diagrams contributing to the /i — )■ 67 decay in the SM (a) and in 
the MSSM (a, b, c). Tilde denotes the superpartner of the corresponding particle. 



Considering only the main contribution stemming from the electromagnetic penguin 
operator the /i — )■ 67 transition amplitude is given by 



r(/i ^ 67) 



327r4 



\Aw± +Ah± +A^±\ 



where the three contributions correspond to the diagrams a), b) and c) in Fig. ^ respec- 
tively. The LFV process in the first two contributions involves neutrinos in the loop. For 
the very small neutrino masses < leV, as it follows from eq.(Il]), one gets a negligibly 
small branching ratio. The same statement is valid for the charged Higgs boson whose 
mass is expected to be a few 100 GeV. Thus, one is actually left with the contribution of 
the sneutrino depicted in diagram c). 

The chargino contribution to A^± is given by 



-1 k- 



V2COS/3 



rriw 



m-± 



\Ual\ 





(9) 



where the matrices U and V are the chargino mixing matrices and is the sneutrino 
mixing matrix. In our convention, the sneutrino flavour eigenstate basis is rotated to 
the sneutrino mass eigenstate basis in the same way, as it is in the neutrino sector: 
z//^ = Vi^. In the limit of a minimal LFV the sneutrino mass matrix in eq. (|9]) is 

represented by the PMNS matrix. The functions fi{x) and f2{x) can be written as 

1 



12(1 - 
1 



x] 



2(1 -x) 



J (2x + — 6x'^ + x^ + In x) 
(1 - 4a; + 3x'^ - 2x'^liax) . 



(10) 
(11) 
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Note that in all the cases one has the product of the PMNS matrix elements UniU*^ which 
being summed over i gives zero due to the unitarity of the PMNS matrix. However, this 
product comes with the function / (^^l^/^'^±^ which depends on the mass of the i-th 

sneutrino. If all sneutrinos are degenerate, this function is universal and the result is zero. 
So the whole contribution crucially depends on the splitting in the sneutrino sector. In the 
MSSM with the universal boundary conditions the splitting is achieved via the non zero 
tau Yukawa coupling of the third generation. Hence, effectively, one has the contribution 
of the third generation 




cos 6'i3 sin 6'i.s sin 6' 



13' 



723- 



12) 



If the angle di^ is zero, the whole contribution vanishes. If, on the contrary, it is big, the 
obtained decay width can contradict the experimental data, as it will be clear later. So 
the value of ^13 becomes crucial. 

Since the muon decays to almost 100% as /i — )■ ez/z/, Ttot = r(yU — ez/z/), and the 
branching ratio Br{^ — )■ 67) can be written as 



i?r(/i — > 67) 



T{fi 67) 
r(/i — !■ ez/z/) 



(13) 



which finally gives, according to eqs. ([6]) and ([8]), 



(14) 



At this point a question can arise whether and to which extent our results would 
be sensitive to the contribution of a hypothetical right-handed sneutrino which is not 
indeed taken into account in eq.([9]). In fact, a right-handed sneutrino does not belong to 
the particle content of the MSSM. Therefore, no mixing between left-handed and right- 
handed sneutrinos can occur and the sneutrino mass matrix is given by 



rri: 



+ |M| cos 2/31 



3x3 







3x3 



03X3 
03X3 



(15) 



However, the MSSM can be extended by introducing right-handed neutrinos and their 
supersymetric partners. Further, with the implementation of a seesaw mechanism of 
type I [2TH2S] the neutrino masses and mixing angles can be generated. In the presence 
of Majorana and Dirac mass terms in the Lagrangian the seesaw mechanism provides 
one light and one heavy neutrino mass eigenstate with masses fn^^.^^^^ = rrvj^/Mn and 
"^fhoav ~ respectively, where higher order terms in 1/Mr are neglected, mu is the 
Dirac mass and ~ 10^^ GeV is the Majorana mass. In this MSSM-seesaw model 
there is a contribution from the heavy neutrino mass eigenstate to the W~ and H~ loops 
and from the heavy sneutrino mass eigenstate to the x~-loop- The contribution from 
the heavy neutrino to the LFV transitions is suppressed by a very small mixing angle 



m 



i-'light 



/Mr and, therefore, is completely negligible. 
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Further, we can estimate the impact of the Majorana mass on the chargino mediated 
LFV process. The Majorana mass enters into the LR, RL and RR blocks of the sneutrino 
mass matrix |26|l27]. Diagonahsing the sneutrino mass matrix (see i.e. Ch. 2.1 in ref. [28] ) 
one obtains corrections to the mass of the hghter sneutrino which are of the order of the 
hght neutrino mass. The heavier sneutrino obtains a mass m^^^^^^ = Mr. Its contribution 
to the chargino mediated process is again suppressed by a very small mixing angle given 
above and, in addition, by the fact that the function /2 vanishes for a big argument, 
/2(x — >■ oo) — )■ 0. Hence, the effect of the introduction of the Majorana masses is negligible 
and there is practically only the contribution from the left-handed sneutrinos to the 
chargino mediated /z — ?• 67 transition. The heavy sneutrinos decouple and the low-energy 
sneutrino mass eigenstates are dominated by the pl components. 

Assuming mSUGRA universality conditions at the unification scale the LL block of 
the sneutrino mass matrix can also be affected by the large mass M/j through the renor- 
malisation group running from the unification scale down to the SUSY scale. This leads 
to a LFV caused by renormalisation group effects. The main idea is that the LL block of 
the sneutrino mass matrix is modified through a radiative correction to the soft slepton 
mass matrix m? . In leading logarithmic approximation the correction leads to off-diagonal 
elements in the sneutrino mass matrix which are typically smaller by a factor of ~ 10~ 
in comparison with the diagonal ones. Collecting all the flavour violation effects into 
the PMNS matrix, this additional flavour violation appears as a correction of the PMNS 
matrix elements of an order of ~ 10^^. Hence, the radiative LFV does not lead to any 
sizeable effects in our case of study. For a comprehensive study of the radiative LFV we 
refer to refs. P^[T7t[25H55] and references therein. 

4 Numerical analysis 

We calculate the /i — )■ 67 branching ratio in the MSSM with minimal LFV and mSUGRA 
universality conditions. The boundary conditions of this so-called constrained MSSM 
(CMSSM) imposed on the multidimensional MSSM parameter space imply that at the 
GUT scale all the sleptons, squarks and Higgs bosons have a common scalar mass mo, all 
the gauginos are unified at the common gaugino mass mi/2, and so all the tri-linear terms 
assume a common tri-linear mass parameter Aq. In addition, at the electroweak scale 
one selects the ratio of Higgs vacuum expectation values tan^ and sgn(/i), where is the 
higgsino mass parameter of the superpotential. So one is left with the five-dimensional 
parameter space (mo, mi/2, ^o? tan /3, sgn{^)). For our numerical analysis we fix the value 
of sgn{n) = 1 and the parameters Aq and tan /3 vary for three different points (mo, mi/2)- 
The dependence on 1^4 1 comes from the RG equations for the running sneutrino masses. 

We plot in Fig. 02]) on the left the relation between the branching ratio Br{fi — )■ 67) and 
Aq for different values of tan/3 for the points (mo,mi/2) = (500,500) GeV, (mo,m-i/2) = 
(1500, 250) GeV and (mo, mi/2) = (500, 900) GeV. The point (mo, mi/2) = (500, 500) GeV 
has been chosen so that it is allowed by other processes, i.e., the branching ratios B — )■ Xgj 
and B — )■ /"*"/", the anomalous magnetic moment of the muon as well as by experimental 
limits obtained by direct searches of the Higgs boson and Dark Matter in the Universe [5^ . 
With the point (mo,mi/2) = (1500, 250) GeV we can analyse a scenario with a heavy 
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sneutrino and a light chargino; it corresponds to the beginning of the focus-point region. 
The choice (mo,mi/2) = (500,900)GeV represents the beginning of the so-called co- 
annihilation region which is the opposite scenario. One should have in mind that values 
for mi/2 < 250 GeV are excluded by direct Higgs searches and B — )■ Xg'y even for small 
tau/S. On the other hand, the process B — >■ n'^fi~ is not compatible with small mo. 
We calculate the numerical values of the PMNS matrix elements with sin^ 613 = 0.016 
obtained by a global analysis of all available neutrino oscillation data [12] . 

Our results show that for the points (mo,mi/2) = (500, 500) GeV and {mo,mi/2) = 
(500, 900) GeV the branching ratio Br{^ — )■ 67) grows with Aq for all values of tan /3 
while in the case of the point {mo,mi/2) = (1500,250) GeV we see the opposite trend. 
The current experimental upper bounds for Br{fi — )■ 67) represented by two horizontal 
solid lines on the plots do not allow tan /3 > 35 even for a small Aq and heavy chargino. 
An increasing Aq together with a small chargino mass imposes even more stringent bound 
on the maximal value of the parameter tan/3. 

Note, however, that these conclusions are valid for the fixed value of sin^ ^13 = 0.016. 
At the same time, as was already demonstrated in eq. f|T2|) . our results crucially depend 
on the neutrino mixing parameter 613. This proportionality is explicitly shown in the plots 
on the right side of Fig. ([3]) where the value of Aq is fixed, Aq = 0, and we treat sin6'i3 
as a free parameter. The branching ratio Br{fi — )• 67) grows up with ^13 and exceeds the 
experimental upper bounds even for small values of tan /3 if sin ^13 is big enough. On the 
contrary, for small values of sin ^13 all the values of tan /3 are allowed. 

5 Discussion 

We have shown that the /i — )■ 67 decay might serve as a direct manifestation of physics 
beyond the SM, in particular supersymmetry. Experimental bounds are very close to 
the predicted values. One has a unique combination of SUSY predictions with the most 
involved measurement in neutrino physics, namely, the mixing between the first and the 
third generations 613. Provided the value of ^13, the yU — )■ 67 decay would place more 
stringent bounds on the MSSM parameter space than other rare decays. One can see 
that the high tan /3 scenario of the MSSM favoured by dark matter abundance |35] might 
contradict the — )■ 67 decay for the essential part of the bulk region if the value of 6*13 is 
big enough. On the contrary, global analysis of the SUSY parameter space surely prefers 
small values of 613. It seems that the resolution of both puzzles might come together. 
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Figure 3: The branching ratio for the /i — > 67 decay as a function of Aq and tan/3 (left) 
and as a function of sin6'i3 and tan/3 (right) for the points {mo,mi/2) = (500,500) GeV, 
(mo,mi/2) = (1500, 250) GeV and (mo,mi/2) = (500, 900) GeV, respectively. From top 
to bottom the curves correspond to tan/3 = 40, 30, 20, 15, 10. The upper solid horizontal 
line corresponds to the experimental upper bound obtained by the MEG experiment 
Br{ii — )• 67) < 2.8 ■ 10"^^ [6] and the lower horizontal line to the one obtained by the 
MEGA experiment Br{fi — )■ 67) < 1.2 ■ 10^^^ [7J. The left dashed vertical line on the 
plots on the right corresponds to the mean value of 6*13 obtained by a global analysis of all 
available neutrino oscillation data sin^ ^13 = 0.016 fT3] and the right dashed vertical line 
represents the upper bound sin^ 6*13 < 0.056 at 99.73% C.L. [2]. For all plots sgni^jj) = 1. 
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